Microquasar LS 5039: a TeV gamma-ray emitter and 
a potential TeV neutrino source 1 



F Aharonian 1 , L Anchordoqui 2 , D Khangulyan 1 and T Montaruli 3,4 

1 Max-Planck- institut fur Kernphysik, Saupfercheckweg 1, 69117 Heidelberg, Germany 
2 Department of Physics, Northeastern University, Boston, MA 02115 
3 Department of Physics, University of Wisconsin, Madison WI 53706 
4 Universita di Bari, Via Amendola 173, 70126 Bari, Italy 

E-mail: f elix . aharonian@mpi-hd . mpg . de 

Abstract. The recent detection of TeV 7-rays from the microquasar LS 5039 by HESS is one 
of the most exciting discoveries of observational gamma-ray astronomy in the very high energy 
regime. This result clearly demonstrates that X-ray binaries with relativistic jets (microquasars) 
are sites of effective acceleration of particles (electrons and/or protons) to multi-TeV energies. 
Whether the 7-rays are of hadronic or leptonic origin is a key issue related to the origin of 
Galactic Cosmic Rays. We discuss different possible scenarios for the production of 7-rays, 
and argue in favor of hadronic origin of TeV photons, especially if they are produced within 
the binary system. If so, the detected 7-rays should be accompanied by a flux of high energy 
neutrinos emerging from the decays of tt^ mesons produced at pp and/or pj interactions. The 
flux of TeV neutrinos, which can be estimated on the basis of the detected TeV 7-ray flux, taking 
into account the internal 77 — » e + e~ absorption, depends significantly on the location of 7-ray 
production region(s). The minimum neutrino flux above 1 TeV is expected to be at the level 
of 10~ 12 cm _2 s _1 ; however, it could be up to a factor of 100 larger. The detectability of the 
signal of mufti- TeV neutrinos significantly depends on the high energy cutoff in the spectrum 
of parent protons; if the spectrum of accelerated protons continues to 1 PeV and beyond, the 
predicted neutrino fluxes can be probed by the planned km 3 -scale neutrino detector. 



1. Introduction 

During the next decade, two km 3 scale underwater/ice neutrino telescopes will start operation, 
both in the Northern (NEMO/km3NeT [Q) and Southern (IceCube j2]) hemispheres. Theoretical 
and phenomenological studies of recent years show that the sensitivities of these detectors should 
allow meaningful probes of different nonthermal galactic and extragalactic source populations 
in high energy neutrinos (see e.g. [Hill])- 

TeV neutrinos are unique carriers of unambiguous information about hadronic processes in 
cosmic accelerators. They are produced in pp or pj interactions through the decay of secondary 
charged pions. Since generally these processes also result in TeV 7-rays of a comparable rate, 

1 The material of this paper was one of the highlight topics of the Plenary talk of F.A. Aharonian on the theory 
of TeV gamma-ray sources at the TAUP05 Conference. Instead of a brief overview of models of all classes of 
detected gamma-ray sources, in this article we decided to focus on high energy processes in the microquasar 
LS5039 which presents a great interest for the Astroparticle Physics community in general, and for gamma- and 
neutrino astronomers, in particular. 



the best candidates to be discovered by planned km 3 scale neutrino telescopes are sources of 
7-rays with fluxes above 1 TeV at the level of (3 — 10) x 1CP 12 ph/ cm 2 s. 

Presently more than three dozens of sources of both galactic and extragalactic origin 
are established as TeV 7-ray emitters. Some of them, in particular four galactic objects 
- RXJ1713. 7-3946 and RXJ0852.0-4622 ( shell type supernova remnants) and Crab Nebula 
and Vela X (plerions), as well as three variable extragalactic (BL Lac) objects - Mrk 421, 
Mrk 501, and 1ES1959+650, do exhibit persistent or episodic fluxes of TeV 7-rays exceeding 
J-y(E > ITeV) = lO^ 11 phcm -2 s~ x . Thus, if a significant fraction of TeV 7-ray fluxes of these 
sources are of hadronic origin, and the spectra of parent protons extend to PeV energies, IceCube 
and NEMO/km3Net should be able to detect the neutrino counterparts of TeV 7-rays. 

Moreover, unlike neutrinos, TeV 7-rays are fragile and suffer strong internal and/or external 
absorption due to interactions with ambient IR/optical photon fields. Therefore, for sources with 
heavily absorbed TeV emission, the chances of detection of counterpart neutrinos could be quite 
high, even for relatively weak 7-ray sources with a TeV flux below 10 -12 ph/cm 2 s. In this regard, 
the microquasar LS 5039 [5] recently detected in TeV 7-rays by the HESS Collaboration jS] is 
of great interest and a promising source also for neutrino astronomy. 

2. High Energy Nonthermal Phenomena in X-ray Binaries 

In one of the first attempts to classify the potential VHE 7-ray emitters, X-ray binaries were 
attributed to the category of serendipitous sources, i.e. objects which "do not have any firm 
a priori basis for their selection as candidate VHE 7-ray sources" []"". Nevertheless, this very 
class of objects played a very important role in gamma-ray astronomy, in particular initiating 
in the 1980s a renewed interest in ground based 7-ray observations. But ironically, the same 
sources ultimately raised questions about the credibility of the results from that controversial 
period of ground-based 7-ray astronomy. As a result, since the early 1990s X-ray binaries have 
not been considered as primary targets for VHE 7-ray observations. This pessimistic view was 
largely supported also by the belief that X-ray binaries should be treated, first of all, as thermal 
sources effectively transforming the gravitational energy of the compact object (a neutron star 
or a black hole) into thermal X-ray emission radiated away by the hot accretion plasma. 

After the discovery of galactic sources with relativistic jets dubbed Microquasars [8], the 
general view on X-ray binaries has dramatically changed; the observations of synchrotron radio 
flares from microquasars demonstrated non- negligible role of particle acceleration (MeV/GeV 
electrons) in these accretion-driven objects. The discovery of 7-rays from LS 5039 provides 
a clear indication that nonthermal phenomena of particle acceleration and radiation in X-ray 
binaries effectively extend into the TeV and, perhaps, also PeV energy regimes. In this regard, 
many original ideas and models (see e.g. [9*1 1101 ITTT 1121 which were initiated by the early 
(wrong ?) claims of TeV and PeV signals from Cyg X-3 and Her X-l (for a review see [7j), 
remain very attractive and applicable, after some modifications, to Microquasars as well. 

Variable TeV radiation has been predicted from microquasars in different astrophysical 
scenarios with involvement of both leptonic [H] and hadronic [H"] interactions. Gamma-rays of 
leptonic [T""j] and hadronic |17| I18j origin associated with termination of jets in the interstellar 
medium are expected also from extended regions surrounding microquasars. 

Hadronic models are attractive in the sense that they predict, from X-ray binaries in general, 
and from microquasars in particular, significant fluxes of both TeV 7-rays and neutrinos due 
to proton-proton, photomeson and photo-disintegration processes that may take place in the 
wind and/or the atmosphere of the normal star, in the accretion disk, and in the jet (see 
e.g. [9t [TT^ [TT1 IT2| ITT""! Hl^ f2?H [2T1 122^ [2"""] |24| j . 

While the models of production of TeV neutrinos demonstrate the feasibility of several 
different scenarios of particle acceleration and interactions in microquasars, one should not 
overestimate the prediction power and robustness of calculations of neutrino fluxes based merely 



on model assumptions such as (i) the total nonthermal energy of the source, (ii) the efficiency 
of the particle acceleration, (iii) the energy spectrum and the maximum energy of accelerated 
protons, (iv) the efficiency of conversion of proton energy into neutrinos. Obviously, the flux 
predictions based on many model assumptions and parameters does not guarantee any reliable 
estimate of neutrino fluxes. While this concerns (perhaps, with some exceptions) all astrophysical 
source populations, in the case of microquasars several model assumptions together result in flux 
uncertainties as large as orders of magnitude. Moreover, in microquasars particles can be boosted 
to multi-TeV energies only if the acceleration proceeds at the theoretically highest possible rates 
allowed by classical electrodynamics. It is not obvious, however, that this should work, from 
first principles, in microquasars. 

The detection of TeV 7-ray emission from LS 5039 has dramatically changed the status of 
model predictions. First of all, the HESS discovery provides the first unambiguous evidence for 
presence of multi-TeV particles in microquasars. Below we demonstrate that if gamma-radiation 
is produced within the binary system, then (2) electrons hardly can be accelerated to multi-TeV 
energies to explain the observed TeV 7-ray signal, and, therefore, (3) the parent particles of 7 
-rays should be protons or nuclei. If so, (4) 7-radiation should be accompanied by TeV neutrinos 
with a production rate comparable to the gamma-ray production rate. Furthermore, since the 
7-rays suffer unavoidable absorption or cascading, (5) the flux of neutrinos should significantly 
exceed the observed flux of 7-rays, and therefore (6) can be (marginally) detectable by the km 3 
class high energy neutrino telescopes. 

3. The case of LS 5039 

While the HESS discovery of TeV 7-ray emission from LS 5039 does support the early theoretical 
speculations of acceleration of particles in microquasars to TeV/PeV energies, the reported data 
do not allow us to relate the 7-ray production region(s) to specific sites of this complex system. 
The upper limit on the source angular size of about 50 arcsec 6 implies that, for a source 
distance d ~ 2.5 kpc [25], 7-ray production takes place within a radius ~ 0.6 pc around LS 5039. 
Although the observed flux does not allow an unequivocal conclusion about the variability of the 
TeV source (Hj, hints have been recently reported [23] of a possible correlation of the TeV flux 
with the 3.9 day orbital period of the binary system. If confirmed by future observations, this 
would be an indication that TeV 7-rays are produced in a compact region, presumably inside the 
binary system. In what follows we will assume that 7-rays are produced within R ~ 10 13 cm, 
which automatically implies very strong 77 absorption. The effect of this absorption should have 
a strong impact on the detectability of TeV neutrinos from this source, provided that detected 
7-rays are of hadronic origin. 

TeV 7-ray production generally requires two components: (a) an effective accelerator of 
particles, electrons and/or protons, up to 10 TeV and beyond; (b) an effective target (converter). 
While the most likely site for particle acceleration in LS 5039 is the jet, which with a speed 
v = 0.2 c and a half-opening angle 6 < 6° extends out to 300 mas (~ 10 16 cm) j2fi| . the bulk 
of 7-rays could be produced both inside and outside the jet. In particular, relativistic electrons 
in the jet can be accelerated through internal shocks [2Z], and consequently high energy 7-rays 
can be effectively produced through inverse Compton scattering [2B1 El El EB EH] ■ However, in 
the jet of LS 5039, because of close location of the very luminous companion star, the effects of 
fast radiative cooling may prevent the electrons from reaching multi-TeV energies. 

The maximum energy of accelerated electrons is achieved when the acceleration time 
approaches the cooling time due to radiative (synchrotron and Compton) losses. The acceleration 
time can be presented in a general form t acc = rjr^/c ~ 0.11 -&rev Bq 1 r\ s, where tl = E/eB is 
the Larmor radius, -E/reV = E/lTeV, and Bq = B/1G is the strength of the ambient magnetic 
field. The parameter 7/ characterizes the efficiency of acceleration: in the case of extreme 
accelerators (maximum possible acceleration rate allowed by classical electrodynamics) rj — > 1, 



whereas for shock acceleration in the Bohm diffusion regime rj ~ 10(v/c)~ 2 (see e.g. [H^)- For a 
star with a temperature kT = 3.5 eV and optical luminosity L± ~ 1.2 x 10 39 erg/s, the energy 
density of the starlight close to the compact object, w r = L+/AttR 2 c, varies between 400 and 
1600 erg/cm 3 at the apastron (R ~ 2.9 x 10 12 cm), and at the periastron (R ~ 1.4 x 10 12 cm), 
respectively |25| . The Compton scattering of TeV electrons in the field of 3kT ~ 10 eV starlight 
takes place in deep Klein-Nishina regime. While in the Thompson regime (E e < 0.1 TeV) the 
cooling time is inverse proportional to the electron energy, ix ~ 0.03i£j^ v s, in the Klein-Nishina 
regime the characteristic Compton cooling time can be approximated with a good accuracy, 
tc ~ Mwq 1 Ejtjy s, where wq = w r /500 erg cm~ 3 . From t acc = tc one finds 

£ c ,max k 2 [B G (v/0.2c) 2 w 1 ] 3 - 3 TeV. (1) 

Formally, for B > 1 G the maximum energy of accelerated electrons can exceed 10 TeV. However, 
for such a large -B-field, the synchrotron losses dominate over the Compton losses. Namely, the 
condition i acc = t sy , where t sy sa AOOBq 2 E^y s is the synchrotron cooling time, gives 

Se, m ax « 3.8 Bq 1/2 (v/0.2c) TeV . (2) 

Equations (0) and (j2J) imply that electrons could be accelerated to energies exceeding 10 TeV 
only in an environment with B < 0.1 G and wq <C 1. Such small fields can be realized beyond 
the binary system, in outer parts of the jet. In this case, 7-rays can be produced through 
the Synchrotron-self-Compton (SSC) scenario |14j . although the contribution from the external 
(starlight) photons could be still significant even at distances quite far from the binary system. 
It was recently argued [HOI that the inverse Compton on starlight cannot explain the hard TeV 
7-ray spectrum, and, therefore, an additional SSC component was proposed to match the HESS 
data. However, the SSC model does not solve the major problem associated with the maximum 
achievable energy of electrons. In fact, the SSC model requires even higher electron energies 
than the external Compton model. 

Thus, any evidence of 7-ray production inside the binary system, e.g. detection of a 
periodic component of TeV radiation would be a strong argument in favor of the hadronic 
origin of these energetic photons. The extension of hard 7-ray spectrum above several TeV 
requires acceleration of protons to > 100 TeV. Within the inner parts of the jet, with a radius 
Rj et ~ 10 7 - 10 8 cm, a magnetic field B > 10 5 G could be sufficient to boost protons up to 
very high energies. The maximum energy is determined by the condition tl < -Rjet> which 
gives E p < 3 x 10 15 (-Rj e t/10 8 cm)(I?/10 5 G) eV. Although these protons are not sufficiently 
energetic to interact with the starlight photons, they can trigger photomeson processes on X- 
ray photons in the accretion disk. These interactions may lead to copious production of high 
energy neutrinos, neutrons, 7-rays and electrons. While neutrons and neutrinos escape the 
source without significant suppression, the electromagnetic fraction of the energy is effectively 
reprocessed, and escape the source mainly in the form of hard X-rays and low energy 7-rays. For 
example, if the luminosity of the accretion disk at UV band exceeds 10 35 erg s _1 , all 7-rays above 
several tens of GeV will be converted into e ± pairs before they can leave the source. Hence, 
the synchrotron radiation of the last generation of electrons will lead to a hard 7-ray spectrum 
oc E , extending to ~ 100 MeV. Interestingly, observations of LS 5039 performed by the 
COMPTEL telescope do show |32j a statistically significant signal of MeV 7-rays with a photon 
index 1.6 ±0.2 and a total energy flux between 1 MeV and 30 MeV of about 9 x 10 -10 erg/cm 2 s. 
If this radiation is indeed related to LS 5039, it would imply a pronounced maximum in the 
spectral energy distribution at the level of 6 x 10 35 erg s . 

Thus, speculating that this radiation is initiated by photomeson interactions (of protons 
with energy > 100 TeV) close to the inner part of the accretion disk, one can estimate the 
expected flux of TeV muon neutrinos at the level of ~ 10 -10 cm -2 s -1 , provided that the 



charged pions decay before interacting with the ambient dense plasma and radiation. In 
the lab frame, the decay time of charged pions responsible for TeV neutrino production is 
t w ± = (E n /m 7T c 2 )T 7r ± « 2.5 x 10^ 3 (E 7T /10 TeV) s. On the other hand, the cooling time of 7r ± 
due to inelastic irp and n'y interactions depends on the ambient gas n p and photon n x densities: 
t-x-p ~ 10 _3 (n p /10 17 cm -3 ) -1 s, and t ny ~ 5 x 10 _3 (nx/10 20 cm _3 ) _1 s. For typical parameters 
characterizing LS 5039, the number density of X-ray photons and the plasma density in the 
region R > 10 7 cm do not exceed 10 20 cm -3 and 10 17 cm~ 3 , respectively. Therefore charged pions 
decay to \i and before interacting with the ambient photons and protons. The production 
of neutrinos from the subsequent muon decay also proceed with high probability as long as the 
magnetic field does not exceed B ~ 10 6 G. This follows directly from the comparison of the 
decay time of muons, t M = (E fM /m fl c 2 )T^ t ~ 0.2(£^/10TeV) s, with their synchrotron cooling 
time t sy « 0.07 ( J B/10 6 G)- 2 (E /t /10 TeV)- 1 s. 

Protons can also interact effectively with the ambient cold plasma, close to the base of the jet 
and/or throughout the entire jet. In what follows we assume that the base of the jet is located 
close to the inner parts of the accretion disk, i.e., the jet axis z is taken normal to the orbital 
plane, with zq ~ 30-Rs (Rs — 3 x 10 5 (MBH/-^©)cm is the Schwarzschild radius). If the magnetic 
field drops as B oc z , the condition of the confinement of protons in the jet, < R, where 
R = Oz is the radius of the jet at a distance z, implies E max oc Sz=constant. Thus, one may 
expect acceleration of protons to the same maximum energy E m&x over the entire jet region. 
However, if there is a faster drop of B with z, the protons at some distance Zt from the compact 
object will start escaping the jet. If this happens within the binary system, i.e. z% < 10 12 cm, 
protons interacting with the dense wind of the optical star will result in additional 7-ray and 
neutrino production outside the jet. 

If the jet power is dominated by the kinetic energy of bulk motion of cold e-p plasma, the 
baryon density of the jet njet can be estimated from the condition Lj et = tt R 2 ct (z) n- }e t(z) m p t> 3 /2. 
The efficiency of 7-ray production in the jet is p 7 = Lj/L p = crppfn J** n- ]e t{z)dz < 1, where L 7 
is the luminosity of VHE 7-rays and L p is the power of accelerated protons. Here, a pp ~ 40 mb 
is the cross-section of inelastic pp interactions, and f n ~ 0.15 is the fraction of the energy of 
the parent proton transfered to a high energy 7-ray photon. Given the recent estimate of the 
black hole mass in LS 5039 M = 3.71}; 3 } M [25]. we set z ~ 3 x 10 7 cm. For the profile of the 
number density we adopt a power law-form, rtj et = no(z/zo)~ s , where s = corresponds to a 
cylindrical geometry, s = 2 to a conical jet, and s ~ 1 is an intermediate case. Expressing the 
acceleration power of protons in terms of the total jet power, L p = nLj ct , one finds the following 
requirement for the jet power, 

L iet « 2 x 10 37 7 ' 34 ' , ' = - erg s" 1 , (3) 

where k is the acceleration efficiency and i> 7i 34 = L 7 /10 34 erg s . The parameter C(s) 
characterizes the geometry /density profile of the jet. For s = 0, 1, 2, one has C(s) = 
zt/zo, ln(z t /zo), and 1, respectively. The case of cylindrical jet provides the highest efficiency 
of 7-ray production. However, since L 7 < l/30Lj ct (assuming w 10% efficiency of proton 
acceleration, and taking into account that the fraction of energy of protons converted to 7-rays 
cannot exceed 30%) the 7-ray production cannot be extended beyond zt ~ 10 4 zo ~ 3 x 10 11 cm. 
The case of a conical jet corresponds to the minimum efficiency of 7-ray production, and thus 
the largest kinetic power of the jet. In this case the bulk of 7-rays are produced not far from the 
base. Finally, in the intermediate case, TeV 7-rays are produced in equal amounts per decade 
of length of the jet, until the jet terminates. 

If 7-rays are indeed produced in pp interactions, one would expect production of high energy 
neutrinos with a rate close to the 7-ray production rate, Q v (E) = £Q 7 (.£?), where £ varies 
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Figure 1. Time averaged 7-ray spectra of LS 5039 formed due to absorption (dotted lines) and 
cascading (solid lines) in the anisotropic radiation field of the normal companion star. Curves 
1 and 2 correspond two different assumptions on the location of the 7-ray production region in 
the jet: at z = 10 8 cm (1), and z = 10 13 cm (2). Dashed curves correspond to the primary 7-ray 
injection rates chosen in a way that the calculated cascade spectra match the observed 7-ray 
flux around 1 TeV. The shaded regions indicate the ranges of MeV/GeV (EGRET) and TeV 
(HESS) 7-ray fluxes 0. 

between 0.5 and 2 depending on the shape of the proton spectrum. However, since 7-rays 
are subject to energy-dependent absorption [S], both the energy spectrum and the absolute 
flux of neutrinos, J v (E) ~ C t J 1 (E) exp[r(E)], could be quite different from the detected 7- 
rays, J 7 (£) « 1.2 x 10~ 12 E^ phcm" 2 s" 1 TeV -1 0. The optical depth t(E) depends 
significantly on the location of the 7-ray production region, and therefore varies with time 
if this region occupies a small volume of the binary system. This may lead to time modulation 
of the energy spectrum and the absolute flux of TeV radiation with the orbital period |33| I34j . 
However, the 77 interactions generally cannot be reduced to a simple effect of absorption. In 
fact, these interactions initiate high energy electron-photon cascades, supported by the inverse 
Compton scattering and 77 pair production processes. The cascades significantly increase the 
transparency of the source. The spectra of 7-rays formed during the cascade development 
significantly differ from the spectrum of 7-rays that suffer only absorption. This is shown in 
Fig. ^ Three processes have been included in the calculations - photon-photon pair production, 
inverse Compton scattering and synchrotron radiation of electrons. The calculations of the 
electromagnetic cascade developed in the photon and magnetic field are based on the method 
similar to the one described in Ref jSS]. Because of the orbital motion, both the absolute density 
and the angular distribution of the thermal radiation of the star relative to the the position of 
the compact object (microquasar) vary with time. In the calculations we take into account 
the effect of the anisotropic (time-dependent) distribution of target photons on the Compton 
scattering and pair-production cross-sections [SHI ■ In Fig. ^ we show the 7-ray spectra averaged 



over the orbital period. 

The solid curves show the cascade spectra, the dotted curves are the spectra of 7-rays which 
the observer would see due to the pure absorption effect (i.e. in the case of effective suppression 
of the cascade), and the dashed curves indicate the injection spectra. Two different locations 
of 7-ray production region have been assumed: (i) close to the base of the jet z < 10 8 cm, and 
(m) well above the base of the jet, z = 10 13 cm. All curves are averaged over the orbital period 
of the system, taking into account the recent data concerning the geometry of the system |25j . 
The cascade spectra are normalized to the reported range of TeV fluxes. They correspond to 
the production rate (luminosity) of 0.1 — 10 TeV 7-rays between 10 (corresponding to the 
dashed curve 2 for d = 2.5 kpc) and 8 x 10 34 erg s _1 (dashed curve 1). Figure ^ shows that the 
cascade 7-rays calculated for z = 10 8 cm agree quite reasonably with the low energy (MeV/GeV) 
data as well, given the large statistical uncertainties and the fact that the EGRET and HESS 
observations are from different epochs. 

Since muon neutrinos should be produced with about the same rate as TeV 7-rays, the 
dashed curves can be used to estimate the neutrino flux expected from the source: J^ M (> 
ITeV) = 1.6 x 10 -11 and 1.9 x 10 _12 cm~ 2 s _1 for the cases (i) and («), respectively. The flux of 
neutrinos could be, in principle, higher if there is suppression of the cascade, as indicated by the 
dotted curves. This could happen if the magnetic field within the binary system exceeds 30 G. 
Although the energy density of this field B 2 /8tt « 35 erg cm 3 is significantly less than the 
energy density of the radiation field, due to the Klein-Nishina effect, the cooling of the electrons 
is dominated by synchrotron radiation in the energy interval E ~ 100 keV — 10 MeV. The flux 
observed by COMPTEL 32 provides an upper bound on the energy released via synchrotron 
radiation, and consequently sets an upper bound on the expected TeV neutrino flux around 
J Ufl (> 1 TeV) = 10~ 10 cm-V 1 . 

Table 1. Number of Vn events with E > 1 TeV expected in 1 yr of observation. We took into 
account that during propagation u^'s will partition themselves equally between z^'s and 17 's 
due to maximal mixing |40j . 
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Future experiments in the Mediterranean Sea, such as ANTARES ^7\, NESTOR [3E|; and 
NEMO can provide meaningful probes of very high energy neutrinos from LS 5039. The most 
advanced project is ANTARES, with construction beginning in the fall 2005 and completion 
expected by early 2007. This detector will have an instrumented area > 0.06 km 2 , with 
an angular resolution better than 1.0° (0.2°) at energies larger than 100 GeV (100 TeV). In 
Table ^ we show event rates expected for ANTARES assuming a power-law neutrino spectrum, 
dN u JdE oc E~ r , with energy cutoff £ max = 10 TeV and 100 TeV, and T =1.5 and 2. We used 
a realistic Monte Carlo simulation to account for the detector response. We also compensated 
for the fraction of the day in which the source is visible through upward-going muon tracks 
by ANTARES, which amounts to about 55% duty cycle (for details seejSHl)- For all four 
combination of parameters T and E max we assume the same energy flux of neutrinos above 
0.1 TeV, .Fe = 10~ 10 erg/cm 2 s. Note that this is comparable with the energy flux corresponding 
to the dashed curve 1 in Fig. ^ and a factor of 5 less than the upper bound on the i^-flux set 
by COMPTEL data. For a 1° search cone, the atmospheric neutrino background is at the level 
of 0.05 yr~ x ; consequently in the most optimistic scenario the full lifetime of the experiment will 
be required to achieve a 5cr discovery. As we also show in Tabled more promising data samples 



will be obtained by NEMO. The larger sensitivity and low background (~ 1 yr _1 ) will allow 
NEMO to probe a broad range of plausible fluxes. 

In summary, our estimates of the neutrino flux associated with recent HESS discovery of 
TeV 7-rays from the direction of the microquasar LS 5039 show that the upcoming neutrino 
experiments will be sensitive to this flux, provided that the observed TeV 7-rays are produced 
in the inner parts of the jet, and that the spectrum of the accelerated parent protons extends 
to PeV energies. 
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